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Abstract 

It has recently been shown (Burer, Math. Program Ser. A 120:479-495, 
2009) that a large class of NP-hard nonconvex quadratic programming prob- 
lems can be modeled as so called completely positive programming problems, 
which are convex but still NP-hard in general. A basic tractable relaxation is 
gotten by doubly nonnegative relaxation, resulting in a doubly nonnegative 
programming. In this paper, we prove that doubly nonnegative relaxation 
for binary quadratic programming (BQP) problem is equivalent to a tighter 
semidifmite relaxation for it. When problem (BQP) reduces to max-cut (MC) 
problem, doubly nonnegative relaxation for it is equivalent to the standard 
semidifmite relaxation. Furthermore, some compared numerical results are 
reported. 
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1. Introduction 



In this paper, we consider the following binary quadratic programming 
problem 

min x T Qx + 2c T x 
(BQP) s.t. afx = bi, z = 1, 2, . . . , m, 

are {-1,1}*, 

where x G R n is the variable, Q G R nxn , c G IT, a { G R n and b { G R 
for all i G / := {1,2, ...,m} are the data. Without loss of generality, Q 
is symmetric, and we assume Q is not positive semidefinite, which implies 
generally that problem (BQP) is nonconvex and NP-hard jlj]. 

Problem (BQP) arises in many applications, such as financial analysis jij, 
molecular conformation problem |3[ and cellular radio channel assignment 
Q. Many combinatorial optimization problems are special cases of problem 
(BQP), such as max-cut problem j5|. For solving this type of problem, a 
systematic survey of the solution methods can be found in Chapter 10 in ^ 
and the references therein. 

It is well-known that semidefinite relaxation (SDR) is a powerful, com- 
putationally efficient approximation technique for a host of very difficult 
optimization problems, for instance, max-cut problem j5|, Boolean quadratic 
program Q. It also has been at the center of some of the very exciting de- 
velopments in the area of signal processing and communications jil, The 
standard SDR for problem (BQP) is as follows: 

min X • Q + 2c T x 
s.t. ajx = bi, Wi G /, 
(SDR) aJXa, = bj, Mi G /, 

Xu = 1, Vz = 1,2, . . . ,n, 

where symbol • denotes the trace for any two conformal matrices. It is 
obviously that problem (SDR) is convex and gives a lower bound for problem 
(BQP) if the feasible set of problem (BQP) is nonempty. Moreover, if the 
optimal solution (x*, X*) for problem (SDR) satisfy X* = x*(x*) T , then we 
can conclude that x* is an optimal solution for problem (BQP). 



Recently, Burer [10J proves that a large class of NP-hard nonconvex 
quadratic programs with a mix of binary and continuous variables can be 
modeled as so called completely positive programs, which are convex but 



2 



still NP-hard in general. In order to solve such convex programs efficiently, 
a computable relaxed problem is obtained by approximation the completely 
positive matrices with doubly nonnegative matrices, resulting in a doubly 
nonnegative programming [11] . which can be efficiently solved by some pop- 
ular packages. For more details and developments of this technique, one may 



refer to [10, lUl, 1121, [13| and the references therein. 

In this paper, a tighter SDR problem and a doubly nonnegative relaxation 
(DNNR) problem for problem (BQP) are established, respectively, according 
to the features of the constraints in problem (BQP) and the techniques of 
DNNP. And, we prove that doubly nonnegative relaxation for problem (BQP) 
is equivalent to the tighter semidefinite relaxation for it. Applying this result 
to max-cut (MC) problem, it is shown that doubly nonnegative relaxation 
for problem (MC) is equivalent to the standard semidefinite relaxation for 
it. Moreover, some compared numerical results are reported to illustrate 
the features of doubly nonnegative relaxation and semidefinite relaxation, 
respectively. 

The paper is organized as follows. In Section [21 a new tighter semidefinite 
relaxation for problem (BQP) is proposed in Section 12.11 Problem (BQP) is 
relaxed to a doubly nonnegative programming problem in Section 12.21 Sec- 
tion [3] and Section 0] are devoted to show the equivalence of two relaxation 
problems for problem (BQP) and problem (MC), respectively. Some conclu- 
sions are given in Section 

2. New relaxation for problem (BQP) 

2.1. New tighter SDR for problem (BQP) 

First, note that problem (BQP) also can be relaxed to the following prob- 
lem by SDR 

min X • Q + 2c T x 
s.t. afx = bi, Wi E I, 
(SDR) afXa, = &?, Vz G /, 

X u = 1, Vz = 1,2, . . . ,n, 

x — xx T y o. 

If the optimal solution (x*, X*) for problem (SDR) satisfy X* = x*(x*) T , it 
holds that x* also is an optimal solution for problem (BQP). On one hand, 
it is worth noting that 

X - xx T h => X h (1) 
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holds always, which further implies that any feasible solution of problem 
(SDR) is also feasible for problem (SDR). It follows that Opt (SDR) < 

Opt (SDR) since the two problems have the same objective functions, where 
Opt(*) denotes the optimal value for problem (*). Therefore, we can con- 
clude that problem (SDR) is a tighter SDR problem for problem (BQP) than 
problem (SDR). 

On the other hand, we can easily verify that the constraint X — xx T y 
is nonconvex, since the quadratic term —xx T is nonconvex. Thus, problem 

(SDR) is nonconvex and not solved by some popular packages for solving 
convex programs. In order to establish the convex representation for problem 

(SDR), a crucial theorem is given below and the details of its proof can be 
seen in Appendix A. 5. 5 Schur complement in (l4j ]. 

Theorem 2.1. Let matrix M G S n is partitioned as 



M 



A B 
B T C 



7/detA ^ ; the matrix H = C — B T A^ 1 B is called the Schur complement 
of A in M. Then, we have the following relations: 
(i) M y if and only if A y and H y 0. 
(it) IfAyO, then M y if and only ifHyO. 



According to Theorem 12.1( h) and p]) , it holds immediately that 



1 X" 

x X 



y o 



X 



xx T y o 



x y o, 



(2) 



i.e., the constraint X— xx T y can be equivalently reformulated as 



1 x T 
x X 



y 



0, which is not only convex, but also computable. So, problem (SDR) is 
equivalently reformulated as follows: 



min X • Q + 2c T x 
s.t. afx = bi, Wi G /, 



(SDR1) 



aJXai = 6?, Vi G /, 



Xu 
1 x 
x X 



1, Wi 

T 



1,2,.. 



n. 



>=0, 
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which is not only convex in form, also can be efficiently solved by some 
popular packages for solving convex programs. 

Here, some examples are given to show that problem (SDR1) is a tighter 
relaxation problem compared to problem (SDR), and the corresponding nu- 
merical results further show that problem (SDR1) is more efficient than prob- 
lem (SDR). These examples are solved by CVX, a package for specifying and 



solving convex programs [15 



Example 2.1. This is a two dimensional nonconvex problem with one linear 
equality constraint, the corresponding coefficients are selected as follows: 



Q 



o 



-3 
-20 



9 



, A= [10, -10], 6 = 0, 



where A = [a x , a 2 , ■ ■ ■ , a m ] T , b = [bi, b 2 , . . . , b m ] T . 

On one hand, we use CVX to solve problem (SDR), then we obtain Opt (SDR) 
— oo ; since problem (SDR) is unbounded below. On the other hand, when 
problem (SDR1) is solved, it follows that Opt(SDRl) = —28 with X = 

in . , . 1]T 



and x = [— 1, 



I I , ... . ... , . Note that the relationship X = xx T holds, 

thus we can conclude that x = [—1, — 1] T also is an optimal solution for 
problem (BQP). The results show that problem (SDR1) is more tighter and 
efficient than problem (SDR) for this problem. 

Example 2.2. This problem is five dimensions with three linear equality 
constraints, the corresponding coefficients are chosen as follows: 







-52 


31 


49 -7 4 








' -20 " 






31 


-16 


-50 -13 -49 






37 


Q = 




49 


-50 


8 44 -30 


c = 




43 






-7 


-13 


44 36 12 






25 






4 


-49 


-30 12 56 






-6 






4 


10 


29 14 -36 " 








11 




A = 




38 


9 


1 -17 23 




b = 




-50 








48 


39 


5 -17 -13 








-36 





If this problem is solved by CVX with problem (SDR), then it returns Opt(SDR) 
— oo since problem (SDR) is unbounded below. When we use problem (SDR1) 
to solve this problem, we have Opt(SDRl) = —307.548, however, the rela- 
tionship X = xx T is not holds for this problem. Therefore, we obtain a 
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tighter lower bound —307.548 for original problem. These results also show 
that problem (SDR1) is more effective than problem (SDR). 

In fact, the constraint x G { — l,+l} n in problem (BQP) further imply 
that the following relationship 

(1 - Xi){\ - xj) > => 1 - Xi - xj + X i:j > 0, VI < i < j < n (3) 

always hold. Combing with ([3]) in problem (SDR1), we get the following new 
semidefmite relaxation problem 

min X • Q + 2c T x 
s.t. ajx = bi, Vi G /, 

afXai = tf, Vi E I, 
(SDR2) X u = 1, Vi = 1,2,..., w, 

1 — ^ — Xj + Xy > 0, VI < z < j < n, 



The above semidefinite relaxation problem (SDR2) is more tighter than 
problem (SDR1) in form, since n ^ n 2 +1 - > inequality constraints are added into 
corresponding problem (SDR2). Furthermore, we will prove that problem 
(SDR2) is equivalent to another convex relaxation problem for problem (BQP) 
in Section [3J 

Now, we test some problems to show that problem (SDR2) is tighter than 
problem (SDR1) from the computational point of view. These problems are 
of one of two types: 



Table 1: Statistics of the test problems 



Type 


Instances 


n 


m 


Function 


RdnBQP 


50 


50 


20 


randn ( • ) 


RdiBQP 


50 


50 


20 


randi ([-10,10],-) 



• RdnBQP. We generate 50 instances of problem (BQP) by MATLAB 
function randn (■)• The symmetric matrix Q is generated by tril (randn (•) , 
-l)+triu (randn (•) ' ,0), and all instances are nonconvex. 

• RdiBQP. 50 instances of problem (BQP) are generated by MATLAB 
function randi ([-10, 10] , •)• The symmetric matrix Q is generated by 



6 



randn(-)+randn(-) '. Each element in the data coefficients is a random in- 
teger number in the range [—10,10]. All instances are nonconvex binary 
quadratic programming problems. 

To compare the performance of two relaxation problems for problem 
(BQP), by using problem (SDR1) and problem (SDR2), respectively, we 
use performance profiles as described in Dolan and More's paper [l6j]. Our 
profiles are based on optimal values for problems (SDR1) and (SDR2). These 
problems are solved by CVX, and the results of performance are shown in Fig- 
ure [TJ From Figure [lj it is obviously that the lower bound which got from 
problem (SDR2) is much greater than that one of from problem (SDR1), for 
test problems RdnBQP and RdiBQP, respectively. Moreover, we find that 
optimal value of problem (SDR2) is strictly greater than that of problem 
(SDR1) for test problems in the experiment. Thus, the performance of prob- 
lem (SDR2) is much better than problem (SDR1) for solving problem (BQP) 
in some sense. 



SDR1 
-SDR2 



SDR1 
-SDR2 



1.01 1.02 1.03 1.04 1.05 1.06 



Figure 1: Left figure is based on optimal values of problems RdnBQP, right figure is based 
on optimal values of problems RdiBQP. 



2.2. Doubly nonnegative relaxation for problem (BQP) 

Recently, Burer [lOj has shown that a large class of NP-hard nonconvex 
quadratic problems with binary constraints can be modeled as so-called com- 
pletely positive programs (CPP), i.e., the minimization of a linear function 
over the convex cone of completely positive matrices subject to linear con- 
straints. Motivated by the ideas, we first establish the CPP representation 
for problem (BQP), and then give its doubly nonnegative relaxation (DNNR) 



7 



formulation. Subsequently, some compared numerical results are presented 
in this section. 

Let z = |(e — x) in problem (BQP), it follows that z G {0, l} n , and then 
problem (BQP) can be equivalently reformulated as follows: 

min Az T Qz - Az T (Qe + c) + e T Qe + 2c T e 

(BQP) s.t. 2aJ z = aje - k, Vi 6 /, 
z G {0, 1}", 

where e denote the vector of ones with appropriate dimension. According to 
Theorem 2.6 in [10] and similar to the analysis in it, problem (BQP) can be 
further equivalently transformed into the following CPP problem 



(CPP) 



min AQ • Z - Az T (Qe + c) + e T Qe + 2c T e 
s.t. 2aJ z = aje — h, Vi G I, 

AafZai = (aje - k) 2 , Vz G /, 
Zu = Z{, vi = 1, 2, . . . , n, 
1 z T 
z Z 



G C 



1+n, 



where C\ +n is defined as follows: 



C 



1+n 



1 G 5 



1+n . 



X 



^z k (z k ) T 



k£K 



U{0}, 



and for some finite {z k } k( z K C i?+ +n \{0}. 

In view of the definition of convex cone in 
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C\ +n is a closed convex 
cone, and is called the completely positive matrices cone. Thus, problem 
(CPP) is a convex problem. However, problem (CPP) is NP-hard, since 
checking whether or not a given matrix belongs to C\ +n is NP-hard, which 
has been shown by Dickinson and Gijen in [17] . Thus, it has to be replaced by 
some computable cones, which can efficiently approximate cone C\ +n . Note 
that the convex cone (S n ) + is self-dual, and so is the convex cone where 
(S n ) + and denotes the cone of n x n nonnegative symmetric matrices 
and the cone of n x n positive semidefinite matrices, respectively. Hence, 
Diananda's decomposition theorem 18] can be reformulated as follows. 



Theorem 2.2. 



C n C S£ n {S n ) + holds for all n. If n < A, then C n 
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By the way, the matrices in H (S n ) + sometimes are called "doubly 
nonnegative" . Of course, in dimension n > 5 there are matrices which are 
doubly nonnegative but not completely positive, the counterexample can be 
seen in 



19 



According to Theorem 12. 2 [ problem (CPP) can be relaxed to the following 
DNNP problem 



(DNNP) 



min AQ • Z - Az T (Qe + c) + e T Qe + 2c T e 
s.t. 2 of z = aje — hi, Wi G I, 

AaJZdi = (aje - ^) 2 , Vi G I, 
Zu — Zi, vi = 1, 2, . . . , n, 
1 ^ T 
z Z 



g 5^ n (Si +n ) + . 



Up to now, the other convex relaxation problem for problem (BQP) is 
established, i.e., problem (DNNP), which is computable by some popular 
packages for solving convex programs, such as CVX, etc. 

Note that problem (DNNP) has n + ra(n 2 +1 ) equality constraints more than 
standard semidefmite relaxation problem (SDR1), and n equality constraints 
more than problem (SDR2), respectively. Thus, the lower bound which get 
from problem (DNNP) is much greater than that one of by problem (SDR1) 
and problem (SDR2), respectively. 

In the following, two types of problems are tested to show the performance 
of problem (SDR1), problem (SDR2) and problem (DNNP), respectively. 
The statistics of the test problems are chosen as follows: 



Table 2: Statistics of the test problems 



Type 


Instances 


n 


m 


Function 


RdBQP 


50 


50 


25 


rand(-) 


RdsBQP 


50 


50 


25 


rands (•) 



• RdBQP. For this type of problems, we generate 50 instances of problem 
(BQP) by using MATLAB function rand(-). The symmetric matrix Q is 
generated by rand(-)+rand(-) and all problems are nonconvex. 

• RdsBQP. The coefficients of 50 instances of problem (BQP) are gen- 
erated by using MATLAB function rands (•), and the symmetric matrix Q 
is generated by rands (-)+rands(-) ' . All instances are nonconvex. 
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We use performance profiles 16( to compare the performance of problem 
(SDR1), problem (SDR2) and problem (DNNP), for problem (BQP), respec- 
tively. The corresponding results of performance are shown in Figure |2j The 
profiles for Figure [2] are based on optimal values of problem (SDR1), problem 
(SDR2) and problem (DNNP), respectively, and these problems are solved by 
CVX. From Figure 121 it is obviously that the performances of problem (SDR2) 
and problem (DNNP) are almost the same, which are better than that one 
of problem (SDR1), for problems RdBQP and RdsBQP, respectively. Thus, 
we can conclude that it is more efficient to use problem (SDR2) and problem 
(DNNP) than problem (SDR1) to solve problem (BQP), from the point of 
view of optimal values. 

Furthermore, we will show the equivalence of the problems (DNNP) and 
(SDR2) in Section [3j 



SDR1 
-SDR2 
-DNNP 



SDR1 
-SDR2 
-DNNP 



Figure 2: Left figure is based on optimal values of problems RdBQP, right figure is based 
on optimal values of problems RdsBQP. 



3. Relationship between relaxation problems 

In this section, we will investigate the relationship between two relaxation 
problems (SDR2) and (DNNP). First of all, the definition of the equivalence 
of two optimization problems is defined as follows. 

Definition 3.1. We call two problems are equivalent if they satisfy the fol- 
lowing two conditions: 

(i) If from a solution of one problem, a solution of the other problem is 
readily found, and vice versa. 

(ii) The two problems have the same optimal value. 
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Now, based on the above Definition 13. 11 the main theorem is given below. 

Theorem 3.1. Suppose that the feasible sets Feas(SDR2) and Feas(DNNP) 
are all nonempty. Then, two problems (SDR2) and (DNNP) are equivalent. 

Proof. The proof can be divided into two parts. First of all, we will prove 
that Opt(SDR2) > Opt (DNNP). 

Suppose that (x*,X*) is an optimal solution of problem (SDR2), let Zy = 
\{l - x* - x* + X*j) and z { = ~(1 - x*), VI < % < j < n, i.e., 

Z = ^(ee T - e(x*) T - x*e T + X*), z = ~(e~ O- (4) 

By ajx* = bi for all i G I and (jl]), we have 

ajx* = aj (e - 2z) = h => 2ajz = aje - bi, Vi 6 J. (5) 
From (]1J) and a[X*ai = bf for all i G I, it follows that 



afX* ai = aJ{AZ - ee T + e{x*) T + x*e T ) ai = bf 
AafZai = (aje -k) 2 , Vz G /. 

Again from (jl]), which imply that 



zz T = \{ee T - e(x*) T - x*e T + X*) - \{e - x*)(e - x*) T 

= \{x* -x*(x*) T ) y o. 



(6) 



Zu = ~(1 - 2x* + X*J = 1(1 - x*) = zt, Vz G /, (7) 

since X* d = 1. 

From 1 - x* - x* + X*. > 0, VI < z < j < n, it holds that 

% > 0,V1 < % < j < n, (8) 
which combining with (JTj), further imply that 

Zi>0,Vl<i<n. (9) 
By Theorem 12.1( h) and (jl|), it follows that 



(10) 



Combining (05) with ©, ©, (|SJ and flgj), it follows that (a, Z) defined by 
01]) is a feasible solution for problem (DNP). 
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Moreover, again from 01]), we have 
AQmZ- Az T (Qe + c) + e T Qe + 2c T e 

= Q • (ee T - e(x*) T - x*e T + X*) - 2(e - x*) T (Qe + c) + e T Qe + 2c T e 
= Q • X* + 2c T x* = Opt(SDR2), 

which further imply that Opt(DNNP) < Opt(SDR2). 

On the other hand, given an optimal solution (z*, Z*) to problem (DNNP), 
and let 

X iS = 1 - 2z* - 2z* + 4Z*., Xi = l- 2z* : VI < i < j < n, (11) 
which imply that 

X u = l-4z* + 4Z^ = l (12) 
since Z* { = z*, Vz = 1, 2, . . . , n. Moreover, 

1-Xi-Xj+ Xij = 1 - (1 - 2z*) - (1 - 2z*) + 1 - 2z* - 2z* + 4Z*. 
= 4Z*. > 0, VI < i < j < n. 

(13) 

From (TTTT) and 2afz* = afe — hi, Vz G /, it follows that 

afx = aj{e-2z*) = b h Vz G /. (14) 
Again from fill I) and 4aJZ*ai = {afe — bi) 2 , Vz G /, we have 



afXcii = af(ee T - 2e(z*) T - 2z*e T + AZ*)a t 
= bl Vz G /. 

From (TTTT) and Theorem 12.1( h). it holds that 

X - xx T = ee T - 2z*e T - 2e(z*) T + AZ* - (e - 2z*)(e - 2z* 
= A(Z* - z*(z*) T ) yo. 



(15) 



(16) 



Bv (TT31. pi]) . (|T51) and (TTBI). we can conclude that (x, X) defined by (JTT]) 
is a feasible solution for problem (SDR2). Furthermore, we have 

X • Q + 2c T a; = (ee T - 2e(z*) T - 2z*e T + AZ*) • Q + 2c T (e - 2z*) 
= 4Z* • Q - 4(z*) T (Qe + c) + e T g e + 2c T e 
= Opt(DNNP), 

which imply that Opt(SDR2) < Opt (DNNP). Summarizing the analysis 
above and according to Definition [3UJ we can conclude that problem (DNNP) 
is equivalent to problem (SDR2). □ 
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Although Opt(SDR2) = Opt(DNNP) in view of Theorem O and Defi- 
nition 13. 1| problem (DNNP) has n equality constraints more than problem 
(SDR2) in form. So, the amount of computation for solving problem (DNNP) 
may be much greater than that one of solving problem (SDR2). In order to 
illustrate this point of view, the compared performance results are shown in 
Figure E] and Figure HI respectively, which are based on the number of itera- 
tions and CPU time for solving problems RdBQP and RdsBQP. The results 
in Figure |3] show that the performance of problem (SDR2) is better than 
that one of problem (DNNP) for problems RdBQP, but the performance of 
problem (DNNP) is better than that one of problem (SDR2) for problems 
RdsBQP, in view of the points of the number of iterations. From the re- 
sults of the performance of CPU time, it is obviously that problem (SDR2) 
is more efficient than problem (DNNP) for solving problems RdBQP and 
RdsBQP, respectively. Summarizing the analysis above, we can efficiently 
solving problem (BQP) by soling problem (SDR2) or problem (DNNP) in 
practice. 




Figure 3: Left figure is based on the number of iterations of problems RdBQP, right figure 
is based on the number of iterations of problems RdsBQP. 



4. An application to max-cut problem 

The max-cut (MC) problem is a kind of important combinatorial opti- 
mization problem on undirected graphs with weights on the edges, and also 



is NP-hard [20[. Given such a graph, (MC) problem consists in finding a par- 
tition of the set of nodes into two parts so as to maximize the total weight 
of edges cut by the partition. 
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Figure 4: Left figure is based on CPU time of problems RdBQP, right figure is based on 
CPU time of problems RdsBQP. 

Let G be an n-node graph, vertex set V := {1, 2, . . . , n}, A(G) the adja- 
cency matrix of graph G, L the Laplacian matrix associated with the graph, 
i.e., L := Diag(y4.(G)e) — A(G). Let the vector u e {+1, — l} n represent any 
cut in the graph G via the interpretation that the sets {i : U{ = +1} and 
{i : Ui = — 1} form a partition of the node set of G, we can get the following 
formulation for (MC) problem 

f\/tn\ max jU T Lu 

s.t. u e{+i,-i} n . 

On one hand, by using the standard semidefinite relaxation technique to 
(MC) problem, we can get the following problem 

max |L • U 

(SDR) s.t. U u = l, Vi = l,2,...,ra, 

uyo. 



Goemans and Williamson [5( have provided estimates for the quality of prob- 
lem (SDR) bound for (MC) problem. By a randomly rounding a solution to 
problem (SDR), they propose a 0.878-approximation algorithm for solving 
problem (MC) based on problem (SDR), which is known to be the best ap- 
proximation ration of polynomial-time algorithm for solving problem (MC). 

On the other hand, according to the technique introduced in Section 
12.21 problem (MC) can also be relaxed to the following doubly nonnegative 
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programming problem 



(DNNP) 



max L • X — x T Le + \e T Le 
s.t. Xa = x { , Vi = 1,2, . . . ,n, 



1 x T 
x X 



e Sf, n n (Si +r 



Remark 4.1. (i) Note that for two relaxation problems (SDR) and (DNNP), 
the feasible sets are all nonempty. It is obviously that the identity matrix E 
is a feasible solution for problem (SDR), and (x,X) = (0,0) feasible for 
problem (DNNP). 

(ii) Compared with problem (SDR), problem (DNNP) has not only n + 
"( n+1 ) new inequality constraints, but also n variables. 

Thus, according to Theorem 13.11 and Remark l4.1( i). we have the following 
theorem. 

Theorem 4.1. Problem (SDR) is equivalent to problem (DNNP). 

Proof. On one hand, suppose that U* is an optimal solution for problem 
(SDR), and let X iS = \{U^ + 1) and x t = |, VI < % < j < n, i.e., 

X = ±(U* + ee T ), x=±e, (17) 

which imply that 

X u = ^(U* i + l)=x l = ^>0,yi<z< n, (18) 
since U^ — l. Then, from U* >z 0, it follows that 



o<^-(^) 2 = i-(^) 2 , 

i.e., 

-1 < U* < 1, VI < i<j <n, 
combining with (JTTl) . we have 



X ij = ^(U* J + l)>0, l<i<j<n. (19) 
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Moreover, from (j!7p . it follows that 

X-xx T = -(U* +ee T ) - ^ee T = hi* h 0, (20) 

it followed by U* h 0. Combining (11%]) . ffP3j) and ([20]) as well as Theorem 
12.1( h). it holds that (x, X) is a feasible solution for problem (DNNP). Again 
from (ED, we have L • X - x T Le + \e T Le = \L • U* = Opt (SDR), which 

further imply that Opt (DNNP) > Opt (SDR). 

On the other hand, suppose that (x*, X*) is an optimal solution for prob- 
lem (DNNP), and let 

U = AX* - 2x*e T - 2e(x*) T + ee T , (21) 

which imply that 

U u = AX*i - 2x* - 2x* + 1 = 1, VI < i < n, (22) 

since X* d — x*, VI < i < n. From (12TI) . it follows that 

U =AX*-2e(x*) T -2x*e T + ee T 

= 4(X* - x*(x*) T ) + (2x* - e)(2x* - e) T h 0. { } 

From (122|) and (|23j) . we can conclude that U defined by (I2T1) is a feasible 
solution for problem (SDP). Furthermore, again from (I2~T|) . we have 

-LmU — L»X* - {x*) T Le + ^e T Le = Opt(DNNP), 

which imply that Opt(SDR) > Opt(DNNP). The proof is completed. □ 

By Theorem 14.11 we can obtain doubly nonnegative relaxation for prob- 
lem (MC) exactly equal to the standard semidefinite relaxation, i.e. problem 
(DNNP) and problem (SDR) are equivalent according to Definition 13. 1[ with- 
out the boundedness assumption of two feasible sets of two problems. 

5. Conclusions 

In this paper, a class of nonconvex binary quadratic programming prob- 
lem is considered, which is NP-hard in general. In order to solve this problem 
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efficiently by some popular packages for solving convex programs, two convex 
representation methods are proposed. One of the methods is semidefinite re- 
laxation, by the structure of the binary constraints of original problem, which 
results in a new tighter semidefinite relaxation problem (SDR2). The other 
method is doubly nonnegative relaxation. The original problem is equiva- 
lently transformed into a convex problem (CPP), which is also NP-hard in 
general. Then, by virtue of the features of constraints in this problem, a 
computable convex problem (DNNP) is obtained through doubly nonnega- 
tive relaxation. Moreover, the two convex relaxation problems are equiv- 
alent. These results are applied to (MC) problem, we can conclude that 
doubly nonnegative relaxation for problem (MC) is equivalent to the stan- 
dard semidefinite relaxation for it. Furthermore, some compared numerical 
results are reported to show the performance of two relaxed problems. 
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